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ABSTRACT. In this paper we establish well-posedness for scalar conser- 
vation laws on closed manifolds M endowed with a constant or a time- 
dependent Riemannian metric for initial values in L°°(M). In particular 
we show the existence and uniqueness of entropy solutions as well as the 
L 1 contraction property and a comparison principle for these solutions. 
Throughout the paper the flux function is allowed to depend on time 
and to have non- vanishing divergence. Furthermore, we derive estimates 
of the total variation of the solution for initial values in BV(M), and we 
give, in the case of a time-independent metric, a simple geometric char- 
acterisation of flux functions that give rise to total variation diminishing 
estimates. 



1. Introduction 

Many physical systems from continuum mechanics can be modelled by 
nonlinear conservation laws. Typically, these partial differential equations 
of hyperbolic type are posed in Euclidean space. For some applications, 
though, the suitable domains turn out to be hypersurfaces or, more generally, 
Riemannian manifolds which, additionally, may change in time. Consider 
for example the shallow water equations on the sphere [30] as a model for 
the global air and water flow. Further examples include the flow of oil 
on a moving water surface, the transport processes on cell surfaces [T| [25]. 
surfactants on the interfacial hypersurface between two phases in multiphase 
flow [TJ, the flow of a fluid in fractured porous media whose fractions are 
considered as a lower dimensional manifold [13], and relativistic flows |13[ 
[26] . Scalar conservation laws have been established as a good model problem 
for studying the nonlinear effects in such systems. 
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In this contribution we consider a nonlinear scalar conservation law posed 
on a closecQ manifold M endowed with a constant or time-dependent Rie- 
mannian metric g on a given time interval /. The conserved quantity 
u : M X I — > R is transported according to a flux function / and un- 
dergoes compression and rarefaction due to the time-dependence of g. Here 
/ = f(x,t,u) is a general family of vector fields on M parametrized by 
time t and by the real parameter u\ in particular we do not assume / to be 
divergence-free. In integral form the conservation law reads 

(1.1) |- / udV+ [ g(f(u),n dK )dV dK = ml 

at JK JdK 

for some smoothly bounded K C M. Here, dV and dVgx are the respective 
volume measures of M and dK, and uqk is the outer unit normal to dK. 
Let the scalar field A : M x / — > R be defined by the relation dtdV = A dV . 
Applying the Gauss theorem to (jl.ip we see that the propagation of the 
conserved quantity u is governed by the Cauchy problem 

d t u + A u + div f(u) = in M x /, 

u(0) = n in M. 

Here, div denotes the Riemannian divergence operator, A accounts for the 
rarefaction and compression due to the time-dependence of the metric, and 
uq : M — > R is a given initial datum. Obviously, this setting includes the 
case of a scalar conservation law posed on a closed, moving hypersurface in 
R rf+1 . This can be seen by pulling back the flux function and the metric 
tensor to the initial hypersurface. We show the existence and uniqueness 
(in the space of measure-valued entropy solutions) of entropy solutions for 
initial values in L°°(M) by the vanishing viscosity method. To this end we 
construct smooth solutions of parabolic regularizations of the conservation 
law and show uniform bounds of the L°° norms of these solutions. Due to 
the time-dependence of the metric and the fact that / is not divergence- 
free the L°° norm might increase with time. With this uniform bound at 
hand, letting the regularization parameter tend to zero we obtain a measure- 
valued entropy solution to the conservation law. The key step in showing 
that the measure- valued solution is, in fact, a uniquely determined scalar- 
valued solution is the proof of an averaged contraction property for arbitrary 
measure- valued solutions, first introduced by DiPerna [8]. The proof we give 
is based on Kruzkov's doubling of variables technique. This approach seems 
to be new in the context of measure- valued solutions and is based on a careful 
application of Lebesgue's differentiation theorem, cf. the discussion in [llj . 
From the averaged contraction property we also deduce the L 1 contraction 
property as well as a comparison principle. In the last part of our paper we 
derive estimates of the total variation of the entropy solution, assuming the 
initial datum to be in BV(M). In doing so it is crucial to work in a global, 
coordinate-independent way in order to get sharp estimates. On a manifold 
the total variation of the solution may increase with time, contrary to the 
situation with constant flux functions in Euclidean space. For a constant- 
in-time metric, though, our estimates allow us to show that a flux function 
gives rise to total variation diminishing estimates if and only if, essentially, 



i.e. compact, without a boundary 
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it is a family of Killing fields. Based on the total variation estimates we 
establish Lipschitz continuity of the solution in time. The validity of total 
variation estimates is important in order to obtain a-priori error estimates 
for finite volume schemes, cf. [181112], Throughout the whole paper we try to 
minimize the regularity assumptions on g and /. In doing so some arguments 
get slightly more involved. To the best knowledge of the authors, the given 
proofs of the strong convergence of vanishing viscosity approximations and 
of the total variation estimates are the first complete ones in the context of 
manifolds, even in the case of a constant-in-time metric. 

The first paper dealing with the existence and uniqueness for nonlin- 
ear transport equations on (time-independent) manifolds was published by 
Panov [23]. In this paper a variant of (|1.2p is considered which is inde- 
pendent of the geometry. Due to this independence Panov can reduce the 
whole problem to the Euclidean case. Recently, numerous beautiful re- 
sults on the theoretical and numerical analysis of conservation on mani- 
folds were published by Ben-Artzi, LeFloch, and their collaborators, e.g. 
[2 El El El [T71 US]. In the seminal contribution [6] a scalar conservation 
law on a time-independent Riemannian manifold with a divergence-free, so 
called geometry-compatible, flux function is considered. The authors present 
results on the existence and uniqueness of entropy solutions, based on the 
averaged contraction property, as well as total variation estimates. There are 
two issues in this paper which need to be clarified, though. The proof of the 
averaged contraction property proceeds as in [51 |2Z] and implicitly assumes 
the existence of very special approximations of general measure-valued so- 
lutions. While the existence of such approximations is a trivial fact when 
dealing with constant flux functions in Euclidean space, this is not so clear 
in the case of manifolds. Furthermore, in the proof of the total variation 
estimates the authors proceed in a local, coordinate-dependent way, and a 
problem occurring then is the fact that the commutator of the Laplacian 
with a non-Killing vector field is a second order differential operator (cf. 
[29], e.g.). Moreover, in [24] Panov proves the existence and uniqueness for 
the Dirichlet problem of a nonlinear transport equation on a manifold with- 
out a Riemannian metric. His approach is based on a kinetic formulation. 
The present paper will appear simultaneously with [10]. In this paper the 
authors show the existence of entropy solutions on moving hypersurfaces 
in M. d+1 (which is a special case of our setting) using a differential calculus 
which is close to numerical analysis. Furthermore, they show the uniqueness 
of these solutions by localizing Kruzkov's doubling of variables technique. 
Finally, they compute approximate solutions using a finite volume scheme. 

The present paper is organized as follows. Notation and preliminary re- 
sults are introduced in Section [21 In Section El we consider the parabolic 
regularization of the conservation law. Section E] then deals with the exis- 
tence and uniqueness of entropy solutions, while in Section El we establish 
the total variation estimates. 
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2. Notation and preliminary results 

For a general introduction to Riemannian geometry see, e.g., [9j[l6]. Let 
M be a smooth, oriented d-dimensional manifold. With a given Riemannian 
metric g = (•,•), assumed to be sufficiently smooth, we associate the Levi- 
Civita connection (or covariant derivative) V = V 9 on the tangential bundle 
TM. In fact, g and V induce a scalar product and a connection on all tensor 
bundles over M. For example, for covector fields a, (3 we havcS (a,p) = 
g 13 otifij, where (g 13 ) denotes the matrix inverse of (gij), and Va(X,Y) = 
X(a(Y)) — a(Vj7) for tangential vectors X,Y. We denote the associated 
norm by \a\ = \a\ g . In coordinates we have ay ; j = (Va)ij = diCtj — 
with the Christoffel symbols Y^-. For a scalar function u, i.e. a 0-tensor 
field, the covariant derivative is simply the differential. In coordinates we 
have u-i = (Vu)i = diU. The vector field gradu is associated with the 
covector field Vu by raising the index, in coordinates uJ = (grad u) 1 = g 13 u-i. 
Furthermore, the adjoint of V is the divergence operator div = div 9 . For 
example, in coordinates for a vector field X we have div X = X 1 ^. By 
V fc we denote the /c-fold application of V which increases the covariant 
order of the tensor field by k. g and V induce Laplace operators acting on 
arbitrary tensor fields. For example, in coordinates for a covector field a we 
have (A g a)k = (Aa)j; = g %3 a^ij = g lJ (V 2 a)jik- We let dV = dV g denote 
the Riemannian volume form, in positively oriented coordinates r E M. d 
we have dV = (det(gij))^ dr. When dealing with objects depending on 
the Riemannian metric (like differential operators, the volume form, etc.) 
we omit the index g whenever this shouldn't lead to confusion. We shall 
make use of Riemannian normal coordinates centred at some point x E 
M. These are given by concatenating the inverse of the exponential map 
exp x = exp 9 with the isomorphism from the tangent space T X M to M. d 
induced by choosing a g-orthonormal basis of T X M. In these coordinates at 
x = we have g^ = 5ij, d^gij = and = 0. For x,y E M we denote 
by d(x,y) = d g (x,y) the Riemannian distance from x to y. We let ^y xy 
be an arc length parametrized, minimizing geodesic connecting x, y which is 
uniquely determined provided that d(x, y) is sufficiently small. Furthermore, 
we denote by P xy : T X M — > T y M the linear isomorphism which is given by 
parallel transport along 7 xy . 

Let U be a coordinate patch of M with smooth coordinate vector fields di. 
Applying the Gram-Schmidt algorithm to these fields we obtain continuous 
vector fields which form a g-orthonormal basis of T X M for each x E U. 
To be more precise, we define inductively 

\di~Ylt=i( d k,ek)e k \ 

These fields obviously have the same regularity as g and are called an or- 
thonormal frame. 



Throughout the paper we use Einstein's summation convention. 
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For a continuous Riemannian metric, all integrable functions /, and al- 
most all x G M we have 

lim-/ \f(y)-f(x)\dV(y) = 0, 

i.e. almost all x G M are Lebesgue points of /. Here B e {x) is the geodesic 
ball of radius e centred at x and j is the mean value integral. This fact can 
be easily proved in local coordinates by using the facts that the Riemannian 
metric is comparable to the Euclidean metric and that in Euclidean space 
for a fixed integrable function almost every point is a Lebesgue point. Now, 
let us state the following simple lemma. 

Lemma 2.1. Let K be a smooth manifold and let M be endowed with a 
continuous Riemannian metric. Furthermore, let f = f(x,k) : M x K — )• 
R be integrable in x for all k and continuous in k for almost all x. We 
assume that for all (xo,ko) there exists a neighbourhood U of xo such that 
the continuity in ko is uniform w.r.t^ all x £ U. Then almost all x G M 
are Lebesgue points for all k G K. 

Proof. We choose a countable, dense subset Kq C K. Obviously, the claim 
is true for all ko G Ko- For arbitrary k G K we have 

/ \f(y, k) - f(x, k)\ dV(y) < i \f(y, k ) - f(x, k )\ dV{y) 

+ { \f(y,k)-f(y,k )\ dV{y) 

JB e (x) 

+ \f(x,k)-f(x,k )\. 

The last two terms on the right hand side get small if we choose ko close 
to k, while for fixed ko the first term tends to zero for almost all x and 
e\0. □ 

We fix a function p G C£°(— 1, 1) such that p > and J Rd p(\x\) dx = 1, 
and for x,y G M we set p t (x,y) := p t (d(x,y)) := \p{ j^r^ )- Let us fix 
coordinates (r l ) centred at x and apply the Gram-Schmidt algorithm to the 
standard basis of R d w.r.t. gij(0). Concatenating the original chart with 
the isomorphism that maps the standard basis to the resulting basis we may 
assume without loss of generality that gij(0) = In these coordinates we 
have 

(2.1) d(0,r) = |r| +G(\r\ 2 ) 

for r — > and g G C 0,1 . In order to prove (|2.1|) let us make a general remark 
that is very easy to prove. If the Riemannian metrics go,gi are comparable 
with constants Co, ci, i.e. Co go < ffi < ci go, then the corresponding distance 

1/2 1/2 

functions are comparable with constants c , c 1 . Now, in the constructed 



with respect to 

In the following we shall employ these coordinates instead of Riemannian normal 
coordinates in order to minimize the regularity assumptions on g. 
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coordinates the metric gij is comparable to the Euclidean metric with con- 
stants of the form l + 0(|r|), depending on the Lipschitz constant of g. This 
proves (|2.ip . 

Furthermore, let us show that in the constructed coordinates the mini- 
mizing geodesic connecting r and points at the origin, up to an error of 
the order 0(|r|), i.e. 

(2-2) 7ro(O) = -n + 0(|r|), 

\r\ 

for r — > and g £ C ' . To this end consider the difference f3 between 7 r o 
and the Euclidean straight line from r to 0, 

P : s i-> 7r0 (s) - (d(0,r) - s)—^—. 

d(0,r) 

As the curve (3 meets the origin at s = and s = d(0, r) we know by the 
mean value theorem that the derivatives of the components f3 l vanish for 
some s l £ [0, d(0, r)]. Taking into account the Lipschitz continuity of /3 we 
can complete the proof of (|2.2p by an application of (|2.ip . 

Let us suppose from now on that the metric depends on time t € I for 
some interval I. If g E C ' 1 , then g(to), g(t±) are comparable with constants 
of the form 1 + O(\to — ti\), and we deduce as above that 

d g (t ){x,y) = d g ( tl )(x,y) + 0{d g(tl) (x,y) \t - h\), 

in particular 

(2.3) d t d g{t) (x,y) = 0(d g{t) (x,y)). 

Of course, the function p e = p e (x,y,t) depends on time, too. Let us again 
fix coordinates centred at x S M. Applying the Gram-Schmidt algorithm 
w.r.t. gij(0,t) as above we obtain time-dependent coordinates such that 
gij(0,t) = 5ij. Using these coordinates we see that for e \ 
(2.4) 

dt / Pe(x,y,t) dV g{t) (y) = d t p e (d(0,r;t)) (det(^(r,i)))5 dr -> 0, 
since dtgij(0,t) = and 

(2.5) d t d(0,r;t) = O(\r\ 2 ) 

for r — > and g G C ' . In order to prove (|2.5p we just note that gij(to), gij(t\) 
are comparable with constants of the form 1 + 0{\r \ |to — *i|)0 

Now, let us define function spaces on M under the assumption that M 
is closed. To this end we fix some finite atlas (ipk,Uk) with a subordinate 
partition of unity (tpk)- We say that a measurable (/,r)-tensor field T € 
X(M) if 7 J;;;;^' w k o <^" 1 G X(<p k (U)) for all k,i n ,j m , where X = L p or 
W k ' p . The norm is defined by 



^It would be slightly easier to prove ()2.4|) with the help of normal coordinates. But 
the technique employed is needed for (|2.3p anyway and, furthermore, the use of normal 
coordinates would need g to be more regular. 
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Obviously, these spaces do not depend on the specific choice of the atlas and 
the partition of unity. Furthermore, the above norms are equivalent to the 
canonical norms w.r.t. the metric g = g(t) (uniformly w.r.t. t on bounded 
time intervals) provided that g is sufficiently regular, e.g. 

imi^ ( M) ~ E f i vir i 2 dv ' 

j=0 jM 

where H k := W k ' 2 . We similarly define the spaces L P (M x I) and L P (M x 
I x M). Whenever we wish to emphasize that given function spaces are 
normed w.r.t. a specific metric we write L p (M,g), etc. 

Finally, we provide a simple technique to produce mollifications of general 
tensor fields. Let T be a measurable (r, Z)-tensor field on M. We define the 
components w.r.t. the chart (ft of the locally supported tensor field T e ^ to 
be 

and we set T e := ^ fc T e \.. Here u e denotes a standard (i-dimensional mollifier 
kernel. T e is a smooth, globally defined (r, Z)-tensor field. By Euclidean 
theory we have T t ->■ T in L P (M) if 1 < p < oo and T £ L P (M), and T e T 
uniformly if T is continuous. Furthermore, for scalar fields T £ L°°(M) we 
have ||T e || ioo / M -j — >■ llTH^oo r M y This follows, like in the Euclidean case, 
from the weak-* lower semicontinuity of the L°° norm and the fact that 
II^IIl°°(M) — II^IIzWmV ^ e similarly define mollifications of tensor fields 
depending additionally on real parameters like time. 



3. Parabolic regularization 

For the rest of the paper we fix a smooth, oriented, closed d-dimensional 
manifold M, a time interval I = (0, T), a time-dependent Riemannian metric 
g on M, and a vector field / on M depending on time and on a real parameter 
u. Throughout this section we assume g, f to be smooth and that / and its 
derivatives of any order are uniformly bounded. In particular, / is globally 
Lipschitz continuous in u. Furthermore, in this section we let the function 
spaces be normed w.r.t. a fixed Riemannian metric, e.g. g(t = 0). We intend 
to construct a smooth, real-valued solution u of the nonlinear parabolic 
equation 

d t u + A u + div f(u) = e An in M x I, 
' 3 ' 1 ' ti(0) = n in M 

for e > and no £ C°°(M). Note that A and the differential operators 
involved depend on the metric and hence on time. It is not hard to check 
that locally A = — g l i dtgij ■ 

Proposition 3.1. For every e > and every uq G C°°(M) (|3.1|) admits a 
unique smooth solution. 
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Proof. We intend to apply the contraction mapping principle. To this end 
we first need to construct appropriate solutions of the linear equation 

dtu + A u + div f(v) = e Au in M x /, 

tt(0) = u in M 

for given v E C(/,L 2 (M)). To this end we set H := L 2 (M), V := H l {M), 
b v :=-divf(v) E L 2 (I,V), and 



a(wo,wi) := / e (dwo , dw\ ) + A wq w\ dV 

J M 

for wq,w\ E V. Note that o(-, •) depends on time since the scalar product, 
A and the volume form do so. Obviously, we have 

\a(w , wi)\ < c \\w \\ v \\wi\\ v , 
(3.2) 2 2 

a(w ,w ) > ci \\w \\ v - A \\w \\ H 

for wo,wi E V and positive constants co,ci,Ao- By Theorem 4.1 and Re- 
mark 4.3 in [T21 Chapter 3] there exists a unique u E L 2 (/,y) admitting a 
distributional derivative u' E £ 2 (/, V) such that u(0) = ujj and 

-<u'(t),w>- +a(u(t),w) =^b v (t),w>- 

for almost all £ E I and all toe!'. Here, -<•,•>- denotes the duality pairing 
between V' and V. For v, v E C(7, the corresponding solutions u, u, and 
almost every t £ I we have 

-< n'(t) - u(t),w y +a(u(t) - u(t),w) =-< -bz(t),w>- . 

Setting w := u(t) — u{t) and applying (|3.2p 9 we get 

~\Ht) ~ u(t)\\ 2 H + Cl\\u(t) - U(t)\\y 

<^b v (t) - ba(t),u(t) - u(t)y +Ao||u(t)-fi(t)||H 

< c \\v{t) - v(t)\\ H \\u(t) - u{t)\\ v + A - u{t)\\ 2 H ■ 

The second inequality follows from the global Lipschitz continuity of /. 
Integrating from to t E / and using the fact that it(0) = u(0) and Young's 
inequality we deduce that 

rt ft 

\\u(t) - u(t)\\ 2 H < c \\v(s) — v(s)\\ 2 H ds + 2\ / \\u(s) - u(s)\\ 2 H ds 
Jo Jo 

< cT\\v- v\\ 2 C (j t H) + 2 j \\u{s) - u(s)\\ 2 H ds. 

By Gronwall's lemma we have 

\\u - u\\ C (i,H) < (°T(1 + 2A re 2AoT ))^ \\v - v\\ c{lH) • 

Replacing / by a smaller interval if necessary we infer that the mapping v i— > 
u from C (I, H) to itself is Lipschitz continuous with a Lipschitz constant 
smaller than 1, hence admitting a unique fixed point. Since the size of the 
time interval is independent of ||«o||# we can repeat this argument a finite 
number of times to obtain a unique solution u on the whole time interval /. 



6 This makes sense due to the embedding {u G L 2 (J, V) \ u G L 2 (I, V')} <-> C{I,H). 
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It remains to show that u is smooth. By a localisation argument this fol- 
lows from Euclidean linear parabolic theory similarly to the proof of Lemma 
2.16 in Chapter 2]. Let us explain the localisation argument by carrying 
out the first step. Let ip G C£°(U) for some coordinate patch U C M. Then 
the localisation ipu solves 

■<(<pu)',w>- + a(<pu, w) 

(3.3) f 

=<b Ul tpw> e / uw Aip + 2w (du,dip) dV 

Ju 

almost everywhere in / and for all w G V. Since u G L 2 (I,V) the func- 
tional defined by the right hand side is in fact induced by some element of 
L 2 (I x U). Thus, writing (|3.3|) in local coordinates we see that ipu solves 
a uniformly parabolic equation in some smooth domain of M. d with vanish- 
ing boundary values, smooth initial datum and right hand side in 1? . By 
Theorem 6.2 in |2Ut Chapter 4] ipu lies in L 2 (I,H 2 (U)) admitting a time 
derivative in L 2 (U x I). As U and <p are arbitrary we see that u lies in 
L 2 (I,H 2 (M)) admitting a time derivative in L 2 (M x /). Using this locali- 
sation argument and differentiating the localized equation repeatedly w.r.t. 
space and time we can adapt the proof of Lemma 2.16 in [2H Chapter 2] to 
show that u is smooth. Note that we need to invoke parabolic LP theory for 
time-dependent operators, see, e.g., [15]. □ 



Adapting a standard argument it is not hard to show that for smooth 
solutions u, u of (13.11) and all < s < t < T we have 



(3.4) / [u(t) - u(t)] + dV g{t) < / [u(s) - u(s)] + dV g{s) 

Jm Jm 

where [v] + := max(u,0) for cel. From this inequality we can deduce an 
L 1 contraction property as well as 

u < u in I x M provided that tt(0) < u(0) in M. 

Note, however, that from this comparison principle we cannot infer uniform 
bounds for the solutions of (|3.ip because constants don't solve this equation. 
This is due to the presence of the factor A and the nonvanishing divergence 
of /. In order to derive uniform bounds we need to adjust the proof of (j3.4|) 
appropriately. 



Proposition 3.2. For every smooth solution u of (J3JJ) and all t G I we 

have 

(3.5) |K*)|| L< » (M) < \\Ml<»{M) eti C2{s)ds + [ efs^l-r c 3 (s) ds, 



o 



where 



'■= — inf X(x,t) — inf d u div f(x, t, u), 

xGM (x,u)eMxM. 



c 3 (t) := ||div/(- ,t,0)|| LOO(M) 
Proof. We set 



for v < 
7] 5 (v) := { v 2 /45 for < v < 25 . 
v — 5 for 25 < v 
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Note that % G W^(M), rf 6 > 0, ^ > 0, 6rfi < 1/2, and for 5 \ and all 

rjs(v) -»• [t;] + , ^(u) u - ^(t;) -> 0, rfg(v) -»• 0. 

We denote the right hand side of (|3.5[) by u. This function solves f)3.1|) with 
u(0) = ||«o|Iloo( M ), A replaced by — C2 and div/(ix) replaced by — C3. Hence 
we have 

d t 775 (ii -u) + r/ s (u -u)(u-u)\ + div (775(1* - u)(f(u) - f(u))) 

- v »(u-u)d(u-u)(f(u)-f(u)) 
= d t r] S (u -u) + T]' s (u - u) {(u - u) A + div(/(n) - /(«))) 
< d t r]s{u — u) + rf s {u -u)(u\ + div/(u) + uc 2 + c 3 ) 
= eAr)j(« — u) — 7/5 (« — u) |d(it — u)\ 2 g < eArjs(u — u). 

Integrating this inequality over (0, t) x M and integrating by parts gives 
f m(u(t) - «(*)) dV g{t) - [ m (u(0) - 2(0)) dV g(0) 

J M J M 

< [ [ U{u-u)d{u-u){f(u)-f{u)) 

J I J M v 

+ A [rf s (u -u)(u-u)- r] S (u - u))j dV g ^ ds. 

Since all integrands are uniformly bounded and converge pointwise by the 
dominated convergence theorem we deduce 

/ [u(t) - u(t)} + dV g(t) < 0. 
J M 

Hence u < u. Replacing ||tto||.r,°o(.M") by — ||wo|Iz,°°(m)' C3 by — C3 and u — u 
by u — u we similarly show that —u<u. □ 



4. Well-posedness 

Throughout this section we assume that g G C 1 ' 1 and f,d u f G C 1 
Furthermore, we assume 5 U div / to be uniformly bounded from below in 
M x I x M. We will establish well-posedness of the scalar conservation law 
(|1.2p for initial values in L°°(M). For arbitrary convex functions r\ : R — > R 
we set 

q(x,t,u) := / r)'(v)d v f(x,t,v) dv 

Jk 

for some feel. For M(R) := (C (R))' we let Prob(R) C M(R) denote the 
subset of probability measures. For a compactly supported v G Prob(R) 
and a continuous function 77 we set (f, 77) := J R r/ dv. 



It might be possible to lower the regularity assumptions. Note, however, that it is 
only for g £ G ' that the Christoffel symbols are Lipschitz continuous and geodesies are 
uniquely defined. Furthermore, ODE theory tells us that for g G C • and x £ M the 
exponential map exp^. is a local C 1 -diffeomorphism. 
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Definition 4.1. For uq G L°°(M) we say that a weak-* measurable^ map- 
ping v : M x / — > Prob(M) C M(R), with supports suppz^f contained in a 
fixed, bounded interval, is a measure-valued entropy solution of (|1.2p i/0 



(4.1) 



d t ip (u, 77) + (p(u,X(ri - t/( ■ ) • ) + divX Q ~ if d[vX f) 

I JM 

+ {v,d<p(q)) dVdt+ [ v?(0)r/M dV g{0) > 
JM 



holds for every convex function rj and every nonnegative test function ip G 
C^(M x [0,T)). We say i/iai u G L°°(M x J) is an entropy solution if the 
family of Dirac measures 5 U is a measure-valued entropy solution. 

The key step in the proof of well-posedness is to show the following av- 
eraged contraction property for measure-valued entropy solutions. In order 
to get clean calculations and a clean result it is crucial to work in a global, 
coordinate-independent way. 

Proposition 4.2. For measure-valued entropy solutions v,a and almost all 
< s < t < T we have 

/ (v t ®o- t ,r})dV g{t) < / {v s ®a s ,r})dV g ( s \, 

JM JM 

where ~n(u, v) := [u — v] + . By interchanging the roles of u and a we see that 
the same inequality holds for r]{u,v) := \u — v\. 

Proof. In (|4.ip we choose rj = r](u, k) := [u — k] + , 

q = q(x, t, u, k) := sgn[u - k} + (f(x, t, u) - f(x, t, k)), 

and 

ip = ip{x, t, y, s) := ip(t) LOi(t - s) p e (x, y, s) 

for fixed k G M, y G M , s G /. Here uji denotes a one-dimensional mollifier 
kernel, p e is the kernel constructed in Section [21 and ip G C^(0, T). We 
integrate the resulting inequality over MxMxI w.r.t. da y)S {k) dV g ^{y) ds. 
Interchanging the arguments of rj, q and the roles of x, t, v and y, s, a we 
obtain an analogous inequality. Adding these two inequalities results in 

ip'(t)lf e (x,t) +ip(t)lf(x,t) + ip(t)lf&t) dV g(t) (x)dt > 0, 

'/ JM 



®More precisely, for each h € L X (M x J,Co(R)) the mapping (x,t) h-> {u x ,t, h(x, t)) is 
measurable. 

®div x denotes the divergence w.r.t. to the explicit dependence on the variable x, i.e. 
fixing the value of u. 
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where 

l{' e (x,t) := / Ui(t- s)p e (x,y,s)(u Xi t®CTy !S ,rj) dV g r s }(y)ds 

J I J M 

I^ e (x,t):=j I u>z{t- s)[d s p e (x,y,s){v Xjt ®a y>s ,rj) 
J I JM v 

+ p e (x,y,s) (u Xtt <g> o- y , s , X(x,t) (r](k 1 ,k 2 ) - dir](k 1 ,k 2 ) h) 

+ A(y,s) (r](k 1 ,k 2 ) - d 2 n{k u k 2 ) k 2 ))) dV glyS) {y)ds 



I e /{x,t):=l I U)z(t- s)(v x . t <3(Ty : t,d x p e (x,y,s)(q(x,t)) 
J I JM 

+ d y p e (x,y,s)(q(y,s)) 

+ Pe{x,y,s) sgn[A;i - k 2 ] + (div f(y, s,h) - div f(x,t,k 2 ))) dV g ^ s) {y)ds. 

Here k±, k 2 denote the integration variables w.r.t. v, a. Let us show that the 
function tends to I\ almost everywhere for e \ 0, where 

If(x,t) : = / p t (x,y,t){u Xtt 0ay t t,v) dV g ( t) (y). 
JM 

To this end we note that 

\ll' e (x, t) — I{(x, t)\ = | / / ojz(t — s) (h(x, s, k±) — h(x,t, k\)) ds du x A 

Jr J i 

< c / -r \h(x, s,k\) — h(x,t,k\)\ ds du Xt t, 
Jr Jt-i 

where 

h{x,s,ki):= I p e (x,y,s) (ay tS ,r)(ki,k 2 )) dV g < s )(y). 

JM 

Since h is continuous in k\ (and x), uniformly w.r.t. s, by Lemma 12.11 the 
mean value integral tends to zero for all k\ and almost all x, t. An application 
of the dominated convergence theorem then proves the claim. We shall use 
this convergence trick several times throughout the proof. Using this trick 
we can similarly show that If tends to 

{v x ,t ® <Tx,t,v) 

almost everywhere for e \ 0. Now, it suffices to prove that I^' 6 and I^' e 
converge to zero almost everywhere if we let e \ and then e \ 0. 

Let us proceed with /|' e (t,x). Noting that r](ki,k 2 ) = d\ri{ki, k 2 ) k\ + 
d 2 r](ki, k 2 ) k 2 and using the Lipschitz continuity of A in space and time we 
can replace the term X(x, t) by X(y, s) to obtain 



Il ,e (x,t) = 0(e + e) + / / u)e(t- s) (v Xjt ®(Ty tS ,r]) 
J I JM 

(d s p e (x,y,s) + p e (x,y,s) X(y,s)) dV g(s) (y)ds. 
By the convergence trick Ifj,' 6 tends to I 2 almost everywhere for e \ 0, where 

JKM) = / Wx,t® °y,t,v) (9tPe(x,y,t) + p e (x,y,t) X(y,t)) dV g r t \(y) 

JM 

+0(e). 
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Since 

(d t p e (x,y,t) + p e {x,y,t) X(y,t)) dV g{t) (y) = d t (p e (x,y,t) dV g{t) (y)) 
and due to (12.311 and (12.411 we have 



\I£(x,t)\= / {{v x .t o-y,t,r)) ~ ( u x,t <8> CTx,t,v)) 

J M 

(d tPe (x, y, t) + Pe (x, y, t) X(y, t)) dV g{t) (y) + 0(e) 
< c -f |(i/ X)t ® 0-^,77) - (i/ X)t (8) (r x ,t,ri)\ dV g (t)(y) + 0{e). 

Another application of the convergence trick shows that the right hand side 
vanishes in the limit e \ 0. 

It remains to show the convergence of I^ e . Using the continuity of div/ 
in space and time we can replace div f(x, t, k 2 ) by div f(y, s, k 2 ), producing 
a term of the order o(l) for e + e \ 0. Furthermore, note that the function 

(.,,^/ ft < w )K„ 



M 



sgn[/ci - k 2 ] + (div/(y,s,fci) - div/(y, s, k 2 ))) dV g ( a )(y) 



is continuous in k\ (and x), uniformly w.r.t. s. Thus we can apply the 
convergence trick to show that tends to ij almost everywhere for e \ 0, 
where 



t) ■= / {v x>t ® o-j, )t , cf> e (x, y, t)(<?(a;, t)) + (Fp e (j;, y, i)(<?(y, t)) 

J M 

+ Pe{x,y,t) sgn[/ci - k 2 } + (div /(y, i, fci) - div /(y, t, k 2 ))) ^ g ( s )(y) + o(l) 
=:I^(x,t)+/| 2 (rE,t) + (l). 

Here I| x and JJ 2 comprise the terms involving <ip e respectively p e . Again 
by the convergence trick we see that I3 2 tends to I3.2 almost everywhere for 
e \ 0, where 

(4.2) 23,2(3:, t) = (v x ,t ® <7x,t, sgn[kt - k 2 ] + (div f(x, t, ki) - div f(x, t, k 2 ))). 

For the rest of the proof we shall omit the variable t. Since grad x d(x, y) = 
— 7i3/(0) an d grad 3 ' d(x, y) = A / xy (d(x,y)), provided that d(x,y) is small 
enough, by the invariance of the Riemannian metric under parallel transport 
we have 

d x Pe{x,y){q{x)) = -d v p e (x,y)(P xy q(x)). 
Furthermore, it's a well known fact from Riemannian geometry that 

f(y, h) ~ Pxyf(x, h) = -^j yx (o)f(y, k) d(x, y) + o{d(x, y)). 

Thus 



Al 



(V(0),V^ (0) (/(y,A:i)-/(y,fe2)))> dV g (y) + o(l). 
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Replacing o~ y with a x produces a term of the order o(l). Indeed, the error 
term may be estimated by 

\(u x ® <r y , h(y,fa,k 2 )) - (v x ® a x ,h(y,ki, k 2 )}\ dV g {y) 

B e (x) 



<cf \(v x ®a y ,h(x,fa,k 2 )) - {v x ® crx, h(x,fa,k 2 )) \ dV g (y) +o(l), 

JB e (x) 

where 

h(y, fa,k 2 ) := sgn[fa - k 2 } + V(/(y, fa) - f(y, k 2 )). 

An application of the convergence trick then proves the claim. Hence, in 
view of (14.21) it suffices to show that 



p'e{d{x,y))d(x,y) (7^(0), vx (o)f{y,h)) dV g (y) -)• div f(x,h) 

M 

for e \ 0. Writing the integral in coordinates centred at x with gij(0) = 5ij 
and taking into account (|2.ip . (|2.2p . and the continuity of V/ results in 

Y,f-A^h) [ ^(M)M^dr + (l). 

ri Jis. d \ r \ 



1 -j 



Using polar coordinates it is easy to check that the integral equals 5 13 . This 
completes the proofF*! □ 

Note that the treatment of I3 1 is the only step in the proof where (so far) 
g and / need to be more regular than merely Lipschitz continuous. 



Lemma 4.3. For a measure-valued entropy solution v we havi 



11 



lim / (u t , I • -u \ ) dV g M = 0. 
t\o Jm 

Proof. We set 

{1 for < t < e 

1- (t- e)/e for e < t < e + e . 
else 

Choosing rj = id and if = f e ,e(t) ip( x ) with an arbitrary ip G C 0,1 (M) in 
(|4.ip and letting e \ we obtain for almost every e £ I 



(u € , id) ip dV g i e ) + / u ip dV g{0) = 0(e). 
m Jm 

By density the left hand side vanishes in the limit e \ for every ip G 
L 1 (M). For r](u) = u 2 we similarly obtain 

(u e , V ) dV g(e) + / rj(u ) dV m > 0(e). 

M JM 



^It would be slightly easier to prove this last step using normal coordinates. But this 
would need g to be more regular. 

^Here and in the following proof one has to take t (resp. e) from the complement of 
some set of vanishing measure. 
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These two facts imply that 



limsup / {i/ e ,r)-T)(uo) - r/(«o)(- - «o)) dV g r e ) < 0. 

e\0 ijU ' v ' 

= 1 ■ -«o| 2 

□ 



Theorem 4.4. For every uq £ L°°(M) f/i.ere exists a unique entropy solu- 
tion u £ L°°{M x I) o/ (|l.'2p . Furthermore, u(t) admits the bound (|3.5p /or 
almost every t £ I . 

Proof. Using the procedure introduced in Section [2] we define mollifications 
u o,9 e ,f € °f u o>9if- Furthermore, we set f*(x,t,u) := f e (x,t,T(u)), where 
the smooth function r : R — > R with |r'| < 1 is the identity for \u\ < N 
and constant for large values of u. Here we choose N to be, say, twice the 
maximum of the right hand side of (|3.5p . From Propositions 13. l l and 13.21 we 
know that there exist smooth solutions u e of (|3.ip admitting the bound (|3.5[) 
with uo,g,f replaced by UQ,g e ,f^. Since this gives a uniform bound there 
exists a subsequence, again denoted by u e , converging weakly-* in L°°(MxI) 
and an associated family v of Young measures o We let 7] : R — > R be a 
smooth, convex function. Note that for the associated entropy flux q e T we 
have 

div 9e q %{tf) = n'(s) div- E d u p T { s ) ds + ??V) du € {d u r T {u € )). 

Multiplying the identity 

d tV (u e ) + r,'(u e ) \ e u e + rf{u e ) div ff e fr(u e ) = eA g * V (u e ) - r? \u e )\du%, 

by a test function ip as in Definition 14.11 integrating over M x I, and inte- 
grating by parts we obtain 

-ff d t <p V (u e ) + <p[\^r,(u*)-r,\ U *)u*) + div^q e T (u e ) 
J i Jm 

- rf(u e ) div^ r T (u')] + M4(V e )) dV g * dt+ [ <p(0) r?K) dV gm 



M 



< 



/ / e A gC ipr](u e ) dV g e dt. 
J i Jm 



Letting e \ we deduce that v is a measure- valued entropy solution of (|1.2|) 
with / replaced by f T . 

Now we let v and a be two measure- valued entropy solutions for the same 
initial datum no £ L°°(M). By Proposition 14.21 for rj(u,v) := \u — v\ and 
almost all < s < t < T we have 

(v t ®o-t,r)) dv g(t) < / ("s,\ • -u \) + {a s ,\ ■ -u \) dV g ^ s) . 
m Jm 



Lemma 14.31 shows that the left hand side vanishes for almost all t £ I, i.e. 
for almost all x, t the measures u x t and a x ,t must concentrate at same point 



12 For the construction of the family of Young measures we may proceed as in the 
Euclidean case, see, e.g., [28] . 
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u(x,t). The resulting function u is an entropy solution which is unique in 
the set of measure- valued entropy solutions. 

By a standard argument from the theory of Young measures we see that 
(u e ) converges to u in L X (M x /). In particular, a subsequence converges in 
L l (M) almost everywhere in /. Hence, in view of the weak-* lower semicon- 
tinuity of the L°°(M) norm and the convergence of the constant H^ollioofAf) 
and the functions C2 = cj 6 , C3 = c| for e \ 0, we deduce that u admits the 
bound (|3.5p with / replaced by f T . Since c\ < C2 and by the choice of N 
the function u, in fact, solves the original equation and admits the original 
bound fl33J. □ 



Remark 4.5. Taking into account the uniqueness of entropy solutions we 
see from the proof that the whole sequence (u e ) of vanishing viscosity ap- 
proximations converges to u strongly in L P {M x I) for all 1 < p < 00 and 
weakly-* in L°°(M x I). Furthermore, for another entropy solution u we 
deduce from Proposition \4-%\ and Lemma \4-3\ the L 1 contraction property 

[ \u(t)-u(t)\ dV g(t) < [ \u(s)-u(s)\ dV g[s) 

for almost all < s < t < T and for s = and almost all t £ I as well as 
the comparison principle 

u < u a.e. in I x M provided u(0) < u(0) a.e. in M. 



5. Total Variation Estimates 

Throughout this section we assume that g,f G Cf oc , and we assume 
d u div / to be uniformly bounded from below in M x / x Furthermore, 
we let the function spaces be normed w.r.t. the time-dependent metric. 

Definition 5.1. The total variation of a function u G L l (M) is defined as 

T Vg(t)(u) ■= sup / u div 9(t) X dV g( t). 
Jm 

The supremum is taken over all smooth vector fields X with \\X\\loo(m g(t)) — 
1. For u & W 1 ' 1 (M) we have 

TV ff(t) («) = / |Vu| fl(t) d^ (t) . 
Jm 

Furthermore, we let 

BV(M) := {u e L l {M) \ TV g(0) (tt) < 00}. 



Remark 5.2. The definition of the space BV(M) does not depend on the 
choice of the metric. To see this, let 31,52 be two metrics of class C 0,1 . Note 



is A 

gain, it might be possible to lower the regularity assumptions (slightly). 
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that there is a unique, positive function v G C 0,1 (M) such that dV gi = v dV g2 . 
Locally we have v 2 = det((gi)ij)/ det((g2)ij)- Hence 



(5.1) 



f u div gi X dV gi = [ u (div g2 (vX) - dv{X)) dV c 
Jm Jm 

+ u (div gi - div g2 )X dV gi . 
Jm 



Since the difference div gi — div 92 is not acting as a derivative on X , the 
modulus of the right hand side of (|5,ip is dominated by 

(5.2) TVg 2 (u) \\v X\\ L oo( Mjg2 ) + c |M| L i( A/)9l ) II^IIl°o(m, S i) • 

This proves the claim. 



CO 
1 — I 

o 

CN 

We shall now derive estimates of the total variation of smooth solutions 
of the parabolic regularization. Again, it is crucial to work in a global, 
coordinate-independent way in order to get clean calculations and sharp 
estimates. 

Proposition 5.3. Assume that g,f,UQ are smooth and let u be the smooth 
solution of (|3.ip for some fixed e > 0. Then we have for all t G I 

TV s(t) (u(t)) < (eXSM* TV gi0) (u ) + f\H *M * c 5 (s) ds) 

exp(e||Ric|| Ll(J)Loo(M)) ), 

where 

c 4 (t) := -inf d t g Xjt (X,X) - inf (V x d u f(x, t, u), X), 

x,X x,X,u 

c 5 (t) :=u max \\V\(- ,t)\\ L i( M) + || sup |Vdiv/(-,t,u)||| Ll(M) . 

«|<«max 

o 

Here u max denotes maximum of the right hand side of (|3.5p and the infima 
are taken over all x G M, all tangent vectors \X\ g ^ < 1, and (for the second 
infimum) all real numbers \u\ < u max . 

Proof. Let us assume for the moment g to be independent of time, in par- 
ticular A = 0. Taking the total covariant derivative of (|3.ip we obtain 
(5.4) 

d t Vu + V div / = e VAu = e (AVu - [A, V]u) = e (AVu - Ric(Vu, • )) . 

Here Ric denotes the Ricci tensor. The commutator identity that is used 
in the last equality follows from the definition of the Riemannian curvature 
tensor R and 

u-ij - u.*ji = u. fc = Ric ifc u. fc . 

In the following we will use the function rjs from the proof of Proposition 
13.21 Taking the scalar product of (|5,4p with |^y^(|Vn|) gives 

d tVs (\Vu\) + (V div f(u),Vu) n 



-- ^ (|Vu|) ((AVu, Vu) - mc(Vu, Vu)) . 



Vu 
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Let us manipulate the second term on the left hand side. We have 

(div/(«)) ;i = = duf^u-j + f, ji {u) + U),J-'i„)u„) :j . 

Concerning the first term on the right hand side of this equation we remark 
that 

duf^u-jui = ((V .ad- 5 u /)(u),grad«) r,' s (\Vu\), 

V 111 I grad u\ 

while with regard to the last term we compute 

(d u f(u)u.i). jU i =div {d u f(u) Vs (\Vu\)) 

+ div(d u f(u)) (\Vu\ r/ 5 (\Vu\) - m(\Vu\)). 

Hence, integrating (|5.5p over M results in 
(5.6) 

^ f r ls (\Vu\)dV + [ ((V sraJ „ d n f)(u),zr a du) n' x (\Vu\) 

(It J y[ J M I grad u 

+ (V gradu div f)(u) rj' s (\Vu\) 

| grad u | 

+ dw(d u f(u)) (\v u \ v' s (\vu\) - m (\v u \)) dv 

= e [ ^fjj^ ((AVu, Vu) - Ric(Vu, Vu)) dV. 

JM |V«| 

Integration by parts (for general tensor fields, see, e.g., [16]) yields 

Jm |Vn| v ' J M \ ' V |V«| )' 

Let us now show that the integrand on the right hand side is nonnegative. 
To this end we let G : R d K denote the convex function x i— y %(|x| K d), 
where |x| Rt j is the standard Euclidean norm. Choosing Riemannian normal 
coordinates centred at some point p G M, the integrand evaluated in p 
equals 

d d 

didjudj(diG){Vu) = ^ didjudkdju{dkdiG)(Vu) > 0. 

ij'=l ij'=l 

Consequently, letting 5 \ in (|5.6p we obtain 
(5.7) 

d f f 

|Vu| dV+ ( ( V Kra d u d»,f) (u) , gr ad u) + (V ffKtM div.f)M 

y | grad ii| | grad u| 



r// 



< -e / Ric (Vu, -— ) 



almost everywhere in / with : = for Vu = 0. An application of 

Gronwall's Lemma completes the proof in the case of a metric independent 
of time. 
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Now, let us drop this assumption on g. Then, instead of (|5.5p we have 
(5.8) 

(*V«, V.) !4flM + <v(. A), V„> + (V div / („), Vu) « 

|Vu| |Vit| |V«| 

^(|Vn|) 



iVd 



((AVu, V-u) - Ric(Vu, Vu)) . 



The last term on the left hand side and the right hand side may be treated 
exactly like in the time-independent setting. The only difference stems from 
the first two terms on the left hand side. But these two terms equal 

d tV5 (\Vu\) + (%)( gradu, q , Srad 7, ) r/ s (\Vu\) + A \Vu\ r,' s (\Vu\) 

zlgradu^ 

+«<VA,^)^(|Vu|) 

= d tV5 (\Vu\) + \ V5 (\Vu\) + (d t g){gmdu, f™ 1 " ) V ' 5 (\Vu\) 

zlgradnl' 

V7; 

+u (VA, — } V ' s (\Vu\) + A (|V«| r/ s (\Vu\) - V s(\Vu\)) . 

The last term tends to pointwise. Hence, integrating (|5.8p over M and 
letting 5 \, 0, instead of (|5.7[) we obtain 

-| / |V«| dF+ /" ((V |S rad„ d„,f)(n),gradn) + (V Kradu div/)(u) 

«<M ./M grad u | grad u] 

+ ^X^"- + " (VA ' M } ^ " " £ I ^ (V "' ^ ^ 

almost everywhere in /. Again, an application of Gronwall's Lemma proves 
the claim. □ 

Concerning the second part of the following theorem, we let v G C 1,1 (M x 
7) denote the unique, positive function such that dV g ( t \ = v(t) dV g ^y 

Theorem 5.4. For every uq G L°°{M )nBV(M) the unique entropy solution 
u G L°°{M x /) of (jl.2p satisfies for almost all t £ I 

(5.9) TV fl(t) (u(t)) < TV 9(0) (n ) + [ c 5 (s)e^^ dT als. 



Furthermore, there exists a representative of u such that for all s,t £ I 

\\u(t)v(t) -u(s)v(s)\\ Ll{Migm < (c 6 + c 7 TV max )|t-s|, 
where TV max denotes the maximum of the right hand side of (|5.9p and 



c 6 := || sup |div/(u)||| Loo(/iL1(M)) , c 7 := || sup |d u /||| Loo(JxAf j. 

\u\<Um 

ax 

\u\<Um 

ax 

In particular, we have u G C 0,1 (J, L X (M)). 

Proof. Using the procedure introduced in Section [2] we define mollifications 
g e ° , f e ° of g,f. One needs to proceed more carefully in order to produce 
mollifications of uq. Let us fix some eo > 0. We let Uq 1 denote the evaluation 
at time ei of the heat flow w.r.t. the metric g € °(0) and the initial value Uq. In 
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particular, we have ||«q 1 |Ii/>° — ll u o||icx>, and the sequence (uq 1 ) tends to mo in 
L X (M). Furthermore, in [22] it is shown that the sequence (TV g e (o)( u o 1 ))ei 
converges to TV 3 e (o)(' u o)- By Proposition 13,11 there exist smooth solutions 
u l ei °f (BHD admitting the bound (|5.3p . with the original data replaced by 
the mollifications. We let first ex, then e, and finally eo tend to zero, the 
distinction between e and eo being necessary in order to avoid a blow-up of 
the Ricci tensor on the right hand side of (|5.3p . So let us start with the 
first limiting process. Inequality (|5.3p shows that the solutions are uniformly 
bounded in L°°(I, W 1 ' 1 (M)), and (|3.ip shows that they admit distributional 
time-derivatives which are uniformly bounded in L°°(I, (W 2 '°°(M))'). Since 

W l ^{M) L 1 (M) (W 2,00 (M)y, 

by the classical Aubin-Lions theorem, see, e.g., [3], a subsequence of (u € ) ei 
converges pointwise almost everywhere to a weak solution of (|3.1|) with 
g, f replaced by g € °,f e °. Taking into account the lower semicontinuity of 
the total variation we easily show that the estimate (|5.3p with g, f replaced 
by g £ ° , / e ° holds for u\ . By another two, very similar applications of the 
Aubin-Lions theorem we let e and eq tend to zero. In doing so we take into 
account the facts that 

TV fl(t) («(t)) < liminf TV g e o(t) (u eo (t)), 
(5.10) e ° 

TV gEO (o)Oo) -> TV s(0) (mo) for e \ 0, 

and the convergence of the functions C4 = c£,c§ = c e 5 ° . Here, the conver- 
gence (|5.10p 9 follows from the argument employed in Remark 15.21 since in 
the estimate analogous to (|5.2p the functions v = v e ° tend to 1 uniformly 
and the constant c = c e ° is of the form O(eo). We thus obtain the first part 
of the theoremQ 

In order to show the second part, let us use in (|4.ip with r/ = id the test 
function 

with and "0 as in the proof of Lemma 14.31 Letting e \ we obtain for 

almost all < s < t < T 

(5.11) 

(u(t) v(t) - u(s) v(s)) i) dV g{0) = / u{t) tp dV g{t) - / u(s) tp dV g{s) 

J M J M 

[ dVdr. 

J M 

Another approximation argument (using again the result from [22]) and the 
chain rule show that the modulus of the right hand side is dominated by 

"* 

\L°°(M) 



(ce \t 



s\ +c 7 J TV 9(t )(m(t)) dr ^ 
Taking in ()5.1ip the supremum over all IIV'IIlo^a/) — 1 completes the proof. 



□ 



14 Of course, instead of invoking the Aubin-Lions theorem we could as well proceed like 
in the last section to deduce strong convergence. 
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Remark 5.5. If we assume that g is independent of time, div/(ti = 0) = 
and 

(5.12) (V x d u f(x,t,u),X)=0 

for all (x,t,u) £ M x I x M and all tangent vectors X £ T X M, then (|5.9[) is 
a total variation diminishing estimate, i.e. 

(5.13) TV g (u(t)) < TV g (u ). 

Note that from (|5.12|) we can deduce that d u div / = duf 1 ^ = 0, in particular 
div / = 0. Condition (|5.12p says that the symmetric part of the differential 
vanishes. Obviously this is equivalent to the claim that the differential is 
skew- symmetric, i.e. 

(V x d u f,Y) = -(V Y d u f,X) 

for all tangent vectors X, Y . Vector fields satisfying this condition are known 
in geometry as Killing fields or infinitesimal isometrics, because the flow 
generated by such a field is a one-parameter group of isometries. 
Obviously, the assumption 

{V x d u f,X) > 

would suffice to obtain (|5.13p . But on closed manifolds this seemingly weaker 
condition is equivalent to (|5.12p . since it implies 

divd u f = (d u fy. t >o. 

By the Stokes theorem this expression must vanish. But then for each fixed 
1 < i < d we have (duf) 1 ^ = 0, which gives (|5.12p . 

Remark 5.6. If we assume g to be time-independent, then, essentially, 
(|5.12p is not only sufficient, but also necessary for f|5. 13|) to hold. Let us 
assume for simplicity that g,f are smooth and that div/ = 0. If (|5.12D 
does not hold, then we have {V xd u f(x,t,u), X) < for some x £ M,t £ 
I,u el, and some tangent vector X £ T X M . We may assume without loss 
of generality that t = 0. Let us fix a positively oriented, orthonormal basis 
of eigenvectors w.r.t. the symmetric part of the differential of d u f(0,x,u), 
with the first eigenvector pointing in the direction of X . We denote by Aj 
the corresponding eigenvalues, in particular Ai < 0. This choice of basis 
induces normal coordinates (r l ) centred at x. We set 

d 

u (r) :=u + KQp(r l /tti) \\p{r l / k 2 ) 

i=2 

for small constants kq > 0, k 2 > n\ > and the bump function p introduced 
in Section® Outside the coordinate patch we set uq = u. For short times 
there exists a unique smooth solution u for these initial values^^ We shall 
choose the constants kq,ki,K2 such that the time- derivative 

^ f \Vu\ dV = - [ ((V grad, d u f)(u),RTadu) dV, 

dt J M J M |gradu| 



'This can be easily shown by the usual implicit function theorem argument. 
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cf. the proof of Proposition 15,31 is positive at t = 0. Note that the right 
hand side is continuous in u and can be written in the constructed normal 
coordinates at time t = as 

(5.14) -f £(A, + O(ko) + 0(k 2 )) T^rP^~ (1 + 0(k 2 )) dr. 

By the change of variables f 1 = r 1 / ' k\ and P = P / 'k 2 for i > 1 we easily see 
that there exists a constant c > such that 

j i f \diuo\ 2 , f \diUo\ 2 j o 

ck,qk 2 < / — - — : — dr < ckq k 2 /Ki and / — - — : — dr<CKoKin 2 

JRd. \VU \] S >d J R d \Vu \ m d 

for i > 1. Thus, taking into account that g(r,0) is comparable to the Eu- 
clidean metric, we deduce that (|5.14j) is positive if 

1/c > (/Co + K2) K 2 /ki + (1 + K + K 2 ) Kl//«2 

/or some new, possibly large constant c > 0. The right hand side can be 
chosen arbitrarily small by setting kq = n 2 = e and Ki = ey^e /or e > 
small enough. 
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